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ABSTRACT 



We present a reconstruction of the mass distribution of galaxy cluster Abell 1689 at z = 0.18 
using detected strong lensing features from deep ACS observations and extensive ground based 
spectroscopy. Earlier analyses have reported up to 32 multiply imaged systems in this cluster, 
of which only 3 were spectroscopically confirmed. In this work, we present a parametric strong 
lensing mass reconstruction using 34 multiply imaged systems of which 24 have newly determined 
spectroscopic redshifts, which is a major step forward in building a robust mass model. In turn, 
the new spectroscopic data allows a more secure identification of multiply imaged systems. The 
resultant mass model enables us to reliably predict the redshifts of additional multiply imaged 
systems for which no spectra are currently available, and to use the location of these systems 
to further constrain the mass model. In particular, we have detected 5 strong galaxy-galaxy 
lensing systems just outside the Einstein ring region, further constraining the mass profile. Using 
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our strong lensing mass model, we predict on larger scale a shear signal which is consistent 
with that inferred from our large scale weak lensing analysis derived using CFh12k wide field 
images. Thanks to a new method for reliably selecting a well defined background lensed galaxy 
population, we resolve the discrepancy found between the NFW concentration parameters derived 
from earlier strong and weak lensing analysis. The derived parameters for the best fit nfw profile 
is found to be C200 = 7.6 ±1.6 and r 2 oo = 2.16 ± 0.10 Mpc (corresponding to a 3d mass equal 
to M200 = 1-32 ± 0.2 x 10 15 hjo M Q ). We find that the projected mass enclosed within the 
Einstein radius for Abell 1689 is M (45") = 1.91 ± 0.27 x 10 14 h 70 M . The large number of new 
constraints incorporated in this work makes Abell 1689 the most reliably reconstructed cluster 
to date. This well calibrated mass model, which we here make publicly available, will enable us 
to exploit Abell 1689 efficiently as a gravitational telescope, as well as to potentially constrain 
cosmology. 

Subject headings: cosmology: observations — galaxies: clusters: individual (A1689) — gravita- 
tional lensing 



1. Introduction 

Galaxy clusters are of considerable cosmological interest, as they are the most recent structures to 
assemble in Cold Dark Matter (cdm) scenarios with the largest angular scale. The detailed distribution 
of mass in galaxy clusters offers a unique opportunity to test structure formation in the Acdm paradigm. 
Cluster mass profiles can be probed via a variety of multi-wavelength data - in the x-ray and optical, SZ 
observations and using measurements of the line of sight velocity dispersion of the cluster member galaxies. 

The lensing effect, i.e. the bending of light by matter along the line of sight from the source to the 
observer, depends only on the mass distribution of the intervening structures, making gravitational lensing 
an ideal tool for measuring the mass profiles of lensing structures. In particular, no additional assumptions 
need to be made with regard to the dynamical state (relaxed or not, in hydrostatic equilibrium or not) or the 
nature (baryonic or not, luminous or dark) of the intervening matter. However, all mass distributions along 
the line of sight contribute to the lensing signal, introducing contamination by foreground or background 
objects. In the core of massive clusters, the surface mass density is well above the critical value enabling 
the use of detected strong lensing features to constrain the inner part of the cluster potential. At larger 
cluster-centric radii, the ellipticities of weakly sheared background galaxies is used to estimate the weak 
lensing effects induced by the cluster potential. Often one limitation of strong lensing is the limited number 
of arcs and multiply imaged systems available to probe the cluster potential. However, with deep Hubble 
Space Telescope (hst) exposure the situation is changing and in this respect Abell 1689 is exceptional since 
it displays the largest number of multiple arc systems ever identified at the center of a cluster. 

In this work, we investigate the mass distribution of the galaxy cluster Abell 1689, using both strong 
and weak lensing. This galaxy cluster is one of the most X-ray luminous z ~ 0.2 Abell cluster, and has the 
largest Einstein radius, around 45", observed to date. Since the size of the Einstein radius is related to the 
size of the critical region of the cluster (i.e. the region where the surface mass density is of the order of the 
critical surface mass density or larger) this cluster was expected to exhibit many multiply imaged systems 
detectable in deep hst exposures. While the field of view of the Wide Field Planetary Camera 2 (wfpc2) was 
too limited to cover the full area of interest, with the Advanced Camera for Survey s (ACS), the situation has 



improved. Indeed, it has provided an unprecedented wealth of arcs in Abell 1689 (jBroadhurst et al.ll2005a . 
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hereafter B05), and many strong lensing studies have been pursued by several groups taking advantage of 
this data-set. Nevertheless, measurements of the mass distribution of Abell 1689 using various techniques 
were not satisfactorily convergent, particularly the strong lensing and weak lensing anal ysis were not found 



quite compatible in terms of the concentration parameter derived from fitting an NFW ([Navarro et al.lll997l ) 
profile to the different data sets. 

We use the hst ACS data in 4 bands: F475W, F625W, F775W and F850LP (see Fig.©. The composite 
color images are crucial to the identification of multiple images, which must have the same colors, see also 
B05 for more details. The new information we use in this work is the spectroscopic confirmation for 24 
multiply imaged systems: we have undertaken a long-term spectroscopic campaign using keck and vlt, 
targeting singly and multiply- imaged sources in Abell 1689. Details of the different observing runs are given 
in Richard et al. (2007), hereafter Paper I. While previous studies claimed 3 spectroscopically confirmed 
multiply imaged systems to constrain their model (and many others with photometric redshifts), we have 
expanded this to 24 spectroscopically confirmed multiply imaged systems. It is the first time a galaxy cluster 
can be modeled using such a wealth of spectroscopic information for arcs. We also use extensive ground 
based imaging: to select cluster members, using data from the ISAAC instrument on the VLT (see Paper I). 
The region covered by ISAAC is a 3' x 3' square centered on the BCG. 

As an additional validation test of our mass model, we will compare the shear profile predicted on 
large scales by the strong lensing model, an d the shear pr ofile measured from CFh12k wide field data. The 
observation and data reduction was done bvlCzoskel hoodi and the photometric and lensing catalog was first 



constructed and described bv lBardeau et al 



Our first aim in this work is to build a reliable mass profile for Abell 1689 and derive key structural 
parameters predicted by such as mass concentration, existence of substructures, halo shape. Secondly, we 
use a new optimization procedure based on the implementation of the Bayesian Monte Carlo Markov Chain 
method in the lenstool software. Thirdly, we present weak lensing measurements from CFh12k data 
that match the predicted shear from our strong lensing model. Hence, thanks to a better identification of 
background sources, we demonstrate reliably that both regimes are found to be in good agreement. Fourthly, 
we report the finding of five strong galaxy-galaxy lenses located at ~ 300 kpc from the cluster center, i.e. 
outside the critical region of the cluster. Such a high rate of strong galaxy-galaxy lensing events clearly 
points to the very high surface mass density of this cluster. 

This paper is organized as follows: first we review in section 2 the results of previous modeling of 
Abell 1689. Section 3 presents our lensing methodology. In section 4, we describe the inclusion of the 
multiply imaged systems. Section 5 reports the new strong galaxy-galaxy lensing systems we discovered. 
Our strong lensing results are presented in section 6 and the weak lensing analysis in section 7. A discussion 
of our results is presented in the concluding section, as well as a discussion on the disagreement found in 
earlier works concerning the concentration parameter of Abell 1689. 

All our results are scaled to the flat, low matter density Acdm cosmology with Hm = 0.3, — 0.7 and 
a Hubble constant Ho = 70 kms -1 Mpc -1 . In such a cosmology, at z = 0.18, 1" corresponds to 3.035 kpc. 
All the figures of the cluster are aligned with WCS coordinates, i.e. North is up, East is left. The reference 
center of our analysis is fixed at the BCG center: RA=13:11:29.52, dec=-01:20:27.59 (J 2000). Magnitudes 
are given in the AB system. 
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Previous Work 



2.1. Optical Spectroscopy and Cluster Dynamics 



The spectroscopic study of the motions of cluster galaxies provides clues to the dynamical state of the 
cluster. It can reveal if the cluster is still und ergoing a merger or if it is already well relaxed (see for instance 



work on the galaxy cluster C10024+1654 bv lCzoske et al 



2001 



2002). 



In early work on Abell 1689. lTeaeue et all |l990h identified 176 cluster members, deriving a high velocity 
dispersion of ~ 2 355 kms -1 . This high value is somewhat misleading, as it reflects the complex dynamical 
state of the cluster. This in fact emphasizes the important r ole that lensi n g can play in the unbiased 
determination of the mass distribution of clusters. Subsequently iGirardi et al.l (|1997l ) . based on 96 redshifts 
in the cluster center, divided the redshift distribution into three distinct subgroups which overlap along 
the line of sight a nd reduced the inferred central velocity dispersion down to ~ 1429 kms' 1 . Recently, 
Lokas et al. (2006!) used s pectroscopic redshif ts to study the kinematics of about 200 galaxies in the cluster 
(mainly coming from the Teague et al.l l|l990h work). They showed that the cluster is probably surrounded 
by a few structures aligned along the line of sight and concluded that the cluster mass cannot be reliably 
estimated only from galaxy kinematics due to the complex structure in Abell 1689. They also find that the 
inferred value of the velocity dispersion depends sensitively on the choice of galaxy sample. Due to technical 
improvements, in particular the advent of high-multiplex spectrographs on 8-10 meter class telescopes, it is 
now possible to get high-qual ity spectra for a large number of cluster galaxies in a short time. The ongoing 
observations bv lCzoskd (|2004l ) using VIMOS on VLT report accurate redshifts for 525 galaxies, spanning from 
the center outwards to 3 h~ x Mpc, providing a deeper insight into the dyna mics of t his cl uster. Using the 
Dressier fc Shectman ( 19881) test to probe for the presence of substructure, Czoskel ( 2004 ) found only one 
apparent distinct group of galaxies that lies ~ 350 kpc to the north-east of the cluster center. The location 
of this substructure corresponds to a group of bright galaxies well identified in optical images of this cluster 
(see Fig. [5]) . The redshift distribution of these galaxies is skewed towards slight ly higher redshifts. On larger 
scales (R > XhT 1 Mpc), no evidence for any substructures is found by the Dressier fc Shectman ( 19881 ) 
test; the outskirts of Abell 1689 look rather homogeneous with no clear signature of any strong clustering. 
The fact that the large scale distribution is regular is also reflected in the velocity dispersion profile which 
decreases from ~ 2100 kms^ 1 in the central bins to ~ 1200 kms -1 in the outer bins (at a cluster-centric 
distance of - 2 hr 1 Mpc). The very high velocity dispersion inferred for the central region can be interpreted 
as a consequence of the non relaxed state of the cluster inner regions, due to an ongoing merger. 

To summarize, optical spectroscopy reveals that Abell 1689 is composed of several structures aligned 
along the line of sight, and that a clear substructure is present in the north-east. 



2.2. X-ray and SZ Studies 



The major assumptions made in most X-ray analysis, in order to derive a mass profile from the observed 
surface brightness profile, is that the intra cluster gas is in hydrostatic equilibrium in the cluster potential. 
This assumption is only appropriate and valid for clusters that have had enough time to relax into equilibrium 
after the last big merger event and that are not experiencing a merger event. 

Abell 1689 has a rather circular X-ray surface brightness distribution, which at firs t look suggests a rela- 
tively relaxe d cluster. However, it ha s no cool core as often observed for relaxed clusters. Andersson 8z Madeiskil 
(j2004h . then I Andersson et all |2006h presented a detailed measurement of the mass from X-ray data obtained 
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with XMM Newton. The temperature map inferred from their analysis shows a clear discrepancy between the 
northern and southern parts of the cluster, with a hint of a temperature gradient in the southwest-northeast 
direction. The redshift map shows a slightly highcr-rcdshift structure to the east at z = 0.185, while the rest 
of the cluster is at z ~ 0.17. These maps imply large-scale relative motion of the intra cluster gas, which 
suggest that the cluster is likely not to be fully relaxed, even if there is no bimodality in the X-ray emission. 
Moreover, a lower than expected gas mass fraction also suggests a complex spatial and dynamical structure. 
To conclude, there is some evidence that Abell 1689 is on the late stage of a merger, but this possible merger 
must be weak given that departures of the gas from a relaxed state are modest. It is com mon to see cluster s 
with a quite regular shape in the x-ray map whereas the lensing mass map looks bimodal j Smith et alj|2005h . 



Furthermore, iBenson et alJ ((200J) measured the Sunyaev-Zeldovich effect in Abell 1689. They found 
that the inferred optical depth of the comptonizing gas might be higher than expected if one considers the 
simple spherically symmetric model obeying hydrostatic equilibrium. They claim that an elongated structure 
along the line of sight would alleviate the discrepancy. 



2.3. Previous Weak Lensing Studies 



Abell 1689 is the first cluster in which a shear signal wa s detected, both in the optical |Tvson et al 

1990) as well as in infra-red observations (jKing et al.l l2002al ). It is also the first cluster for which the 



lens magnificati on effect (depletion curve) has been used to measure the absolute mass ( Taylor et al. 19981 : 
Dve et al.lliool . 



Tyson et al.l (|1990r ) detected the weak lensing effect on a relatively small field of view (720x720 kpc 2 



and selected the background galaxy population by their blue B-R color. This early work concl uded that the in- 
ferred dark matter distribution correlates well with the projected light distribution. Later on. lTvson fe Fischer 
(|1995l ) extended the weak shear analysis to 1 hr 1 Mpc from the cluster center. They found that the surface 
mass density follows a steeper than isothermal profile on scales from 200 hr 1 kpc to 1 h^ 1 Mpc. Comparing 
with N-body simulations in various cosmogonies, they found evidence for the A > CDM model. 



Wide Field ob servations of Abell 1689 were ac quired on the eso/m pg Imager ([Clowe fc Schneiderll2001 
King et al.ll2002bl ) and CFh12k camera at CFHT ([Bardeau et al.l 120051 ). allowing to probe the cluster from 
the center to beyond the virial radius. These three studies used only R band imaging, selecting a background 
galaxy catalog with a magnitude cut-off. The inferred shear profiles appear flat in the inner region, likely 
due to contamination of their background catalogs by cluster members. The Einstein ra dius was estim ated 
to be of the order o f 22" (much smaller than the currently-accepted value of 45"). Clowe (|2003h and 



Bardeau et al.l (|2007l ) later included B and I imaging from CFHT to better select a background galaxy catalog 
with a color cut-off. The derived shear profile is then found to be steeper in the inner part, indicating a 
weaker contamination by cluster members. The M200 mass is found to change l i ttle wit h the color selection, 
but the nfw concentration parameter increases (see Table H]). Broadhurst et al. ( 2005bh from SUBARU multi- 
color data, using a color cut-off to select background members, found a much steeper shear profile than the 
one measured in all previous studies, and a combined fit to this weak lensing data and the strong lensing 
data gives an nfw concentration parameter of 1 0.8, which is not f ully c ompatible with the concentration 
parameter derived from the strong lensing studies. iMedezinski et alJ ([20061 ). using similar SUBARU data, when 
fitting an NFW profil e to the weak lensing data only, finds a concentration parameter as high as 22. Moreover, 
Umetsu et al. ( 2007h used the same SUBARU multi-color data to infer a 2d mass reconstruction from combined 
distortion and magnification data of their background galaxy sample. This study confirmed the Id analysis 
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by Broadhurst et al. ( 2005bh . in particular the high concentration parameter. Recently, Leonard et al. ( 2007 ) 
measured the shape of the background galaxies present in the ACS field of view to estimate the reduced shear 
signal at a radius of 100" from the center of the cluster. They find q = 0.2 ±0.03, a value wh i ch is s maller 
than the one derived from SUBARU data by Broadhurst et al. ( 2005b ) and Medezinski et al. ( 20061 ) (who 
found g ~ 0.4). Note that the contamination by foreground members in this study is likely to be limited, 
since the authors have a much better control of this issue than the ground based studies: they analyze a 
small field of view where most of the cluster members have spectroscopic measurements. Moreover they are 
going much deeper than ground based studies, so the relative contamination is likely to be weaker. This 
study identifies a second mass concentration located somewhere in the north-east part of the ACS field, but 
they were not able to accurately constrain its location. 



2.4. Previous Strong Lensing Studies 



The first strong lensing modeling study of Abell 1689 was done by [Miralda-Escude fc Babul (1995), 
using the detection of giant arcs from ground based data by Tyson et al. ( 1990h to construct a simple mass 
model. This early study already required two mass clumps, a central one and a second one associated with 
the galaxy group in the north-east. They pointed out a discrepancy, of a factor greater than 2, between the 
lensing and x-ray mass estimates, that is likely to arise from the existence of these two clumps. 

The new ACS instrument has allowed strong lensing studies on Abell 1689 to be pursued with greater 
accuracy. The first published strong lensing study of ACS observations of Abell 1689 is the work by B05. 
This remarkable analysis presented the discovery and identification of 30 multiply imaged sources, spanning 
redshifts from z ~ 1 up to 5.5. Subsequent strong lensing works on Abell 1689 have benefited from this 
early identification, although some authors disagree on the identification of a few multiply imaged systems. 
The modeling technique employed by B05 is partly non parametric, and places the initial large-scale mass 
clump aligned with the ligh t distribution, but finds a final mass model which is distributed rather differently 
from the light distribution. [Diego et al.1 (|2005l ) used a fully non parametric method to reconstruct the mass 
using the same data set, where they do not assume that the initial mass distribution follows the light. Their 
results are in agreement with those from B05. 



Zekser et al.l ([20061 ) presented a model where they used the same set of multiple images as B05. In 
this analysis, profiles are used to define the initial form of the galaxy component in agreement with the 
light, but the halo component is allowed to vary using the nfw profile and perturbations added to it. 
These perturbations are modeled as shapelets ([Refregierl 120031 ). adding flexibility in the modeling and thus 
accommodating a broad range of density profiles from cuspy, using a nfw profile to various modifications 
of it. Shapelets are used to add perturbations to the deflecti on field (and not th e mass, which make things 
more complex). Interestingly, the resulting description of the Zekser et al. ( 2006[ ) analysis is dominated by a 
shapelet, representing ~ 85% of the halo mass, whereas the nfw represent only ~ 15%. Thus the shapelet 
contribution to the modeling is more than a perturbation since it dominates the halo budget, suggesting that 
the central part of Abell 1689 is not well described by a nfw profile. They further report some problems 
with a few multiply imaged systems identifications, and they construct models using subsets of multiple 
images, which considerably improves the goodness of the fit (reducing the mean source plane scatter from 
1.19" to 0.74"). 



Recently. lHalkola et al.l (|2006l ). hereafter H06, presented the first fully parametric mass reconstruction 
using the ACS data. They use two large scale dark matter clumps, one associated with the center of the 
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cluster, and the other with the main north-eastern substructure. They report the misidentification for a few 
multiply imaged systems, and propose new ones. One of the key featur es of this study i s the careful modeling 
of the galaxy member perturbation to the main potentials. They use a Brainerd et al. ( 19961 ) profile (bbs) to 
model cluster galaxies, parameterized by a central internal velocity dispersion and a half mass radius. They 
carry out an analysis using the observed fundamental plane in Abell 1689 to derive velocity dispersions for 
176 galaxies with AB magnitudes brighter than 22, and use two different scaling laws to relate the half mass 
radius to the velocity dispersion. They obtain a strong lensing model that is slightly superior to that of B05 
(spatial resolution of RMS of 2.7" compared to 3.2"), and claim the difference is most likely a result of the 
careful inclusion of the cluster galaxies. Moreover, it is worth noting that the RMS of H06 can be regarded as 
superior to previous studies because of the small number of free parameters used in the modeling compared 
to non-parametric approaches. 

In o rder to probe the cluster potential on all scales, both B05 and H06 used weak lensing data from 
SUBARU ( Broadhurst et al. 2005bh in their analysis. When fitting the SUBARU shear profile only with an 
nfw profile, H06 finds a concentration parameter of 30. According to their analysis, for the NFW profile, the 
parameters obtained from strong lensing and weak lensing disagree at the ~ 3cr level. The high concentration 
of an nfw profile fit to weak lensing data is incompatible with both the strong lensing results presented by 
H06 and B05. 

From the above descriptions, it appears that this massive cluster is not fully relaxed. Optical spec- 
troscopy reveals several structures elongated along the line of sight and a sub group of galaxies that lies ~ 
350 kpc to the north-east of the cluster center. Parametric strong lensing studies al l need a second mass 
clump around this galaxy sub-group in their modeling, and the recent flexion study by Leonard et al.l (|2007h 
confirms this result. The absence of cool core also supports a non fully relaxed state for Abell 1689, though 
the X-ray analysis does not support any bimodality as observed in more extreme cases as the bullet cluster 
( Andersson et al. 2006h . Moreover, the large strong lensing cross section observed for Abell 1689, with an 



Einstein radius of about 45 ", also points to some merging processes acting in this cluste r that boost the 
strong lensing cross section (jTorri et alJl2004t iFedeli et alJl2006t iFedeli fc Bartelmannll2007l ) . This likely non 



relaxed state is observed in 70% of clusters at z ~ 0.2 (ISmith et al 



2005) and is not at odd with simulations. 



The very unsatisfactory point is that so far no coherent model was found combining strong and weak lens- 
ing, and more precisely the concentration parameter derived from SUBARU weak lensing data was so high 
compared to the one derived from strong lensing and to expectations of massive clusters that form in N-body 
simulations. Therefore we decided to revisit the lensing modeling of Abell 1689, because such a cluster with 
so many constraints is the best bet to provide a complete coherence between all the modeling efforts of its 
mass distribution. 



3. Lensing Methodology 



We describe below the strong lensing m ass reconstruction that we have used for Abell 1689. Schemat- 
ically, we have use the parametric method ( Kneib et al. 19961) as implemented in the lenstool software 
which is publicly available^. In short we propose that the projected mass is made of two main dark matter 
clumps that are very close to, or possibly in the process of, merging. Then, we have introduced the de- 
flecting mass of individual galaxies using common scaling laws that statistically links the parameters of the 



1 http:// www. oamp.fr/ cosmology /lenstool/ 
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galaxies to their luminosity. We use the observational constraints (multiply imaged systems) to optimize 
the parameters used to describe the mass distribution: this is what we refer to as optimization procedure in 
the following sections. An alternative method is the non-parametric approach that have been used by other 
workers. In Appendix we discuss the main differences between parametric and non-parametric methods. 



3.1. Modeling the Dark Components 



We model the dark matter component separately on small and large scales. The large scale model 
considers clumps that correspond to the more extended cluster mass component, in contrast to smaller-scale 
clumps that are associated with individual cluster galaxies. 

In this work we have modeled all dark matter halo s present in the cluster by a truncated Pseudo Isother- 
mal Elliptical Mass Distribution (piemd, derived from lKassiola fc Kovnerlll993f ) scaled to their appropriate 
m asses: from galaxy hal os to dominant cluster halos. Detailed properties of this mass profile can be found 



Limousin et alj (|2005l ). The density distribution is given by: 



p(r) 



Po 



(l+rVr c 2 orc )(l + rVr c 2 ut ) 



(1) 



This mass profile is parameterized by a central density, p Q , linked to the central velocity dispersion, <j , which 
is related to the depth of the potential well. It is described using two characteristic radii that define changes 
in the slope of the density profile. In the inner region, the profile is described by a core with central density 
po- The transition region (r co re < r < r cut ) is isothermal, with p ~ r~ 2 . In the outer parts, the den sity 

„-2 



progressively falls from p : 



to p : 



introducing a cut-off. Recent work bv lKoopmans et al.1 1)20061 ) has 



shown that early type galaxies are found to be isothermal in their inner parts, with no significant evolution 
with redshift up to ~ 1, supporting the use of an isot hermal model to describe early type galaxies. The PIEMD 
has successfully bee n used to model galaxy clusters ( Kneib et alJll99d: Smith et al. 2005 : Coyone et al. 20061 



2n02ab 



Limousin et al 



Rich ard et al.ll2007t) as well as early type elliptical galaxies (jNataraian et al.lll99S , 
20071 ). Note that when using this profile to model galaxy scale dark matter halos, we fix r corc to be arbitrary 
equal to 0.15kpc, making it very similar to the BBS profile. Then r core is scaled with the luminosity as 
explained below. 

Each clump, of either class, can be fully characterized using seven parameters: the center position (RA, 
dec), the ellipticity e, position angle 8, and the parameters of the density profile, <7o> 



and 



1"cut- 



We consider two large scale dark matter clumps for Abell 1689. The first clump is associated with the 
center of the cluster, and the second clump is associated with the north-east substructure. As a first guess, 
the position is set to coincide with the luminous barycenter of each structure, but during the modeling the 
tolerance limits on the position parameters are large enough to allow the center of the clumps to differ from 
the luminous component. For each clump, we fixed the value of the r cut parameter to reasonable values 
(Table [3]): this parameter describes the properties of the mass distribution on large scales, much larger than 
the radius over which the strong lensing constraints can be found. Thus in practice strong lensing cannot 
give any reliable constraints on this parameter. For the second north-east clump, we decided to limit its 
velocity dispersion to ~ 500 km s^ 1 since no X-ray emission is detected from this region. The total number 
of free parameters describing the large scale dark matter clumps is equal to 12. 

The parameters describing the dark matter halos associated with three individual galaxies were found 
to play an essential role for reproducing multiply imaged systems (systems 1, 2 and 6) because they perturb 
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some arclets (see e. a. iMeneghetti et al.ll2007f ). thus they are also allowed to slightly vary in the optimization 
procedure. The positions will be chosen within ~1 kpc from the visible center of the associated galaxy, and 
other parameters are allowed to vary between reasonable limits during the optimization procedure. Then 
we include the other cluster galaxies in the optimization procedure by applying empirical galaxy scaling 
relations based on the galaxy luminosity in order to reduce the number of free parameters. Most of the 
cluster members are elliptical, but we also consider a few late-type galaxies for which we have measured a 
spectroscopic redshift equal to the redshift of the cluster. Note that we only consider the galaxies that are 
observed inside the hst ACS field of view. 

To extract early-type cluster galaxies, we plot the characteristic cluster red-sequences (R-k) and (b-r) in 
two color-magnitude diagrams and select the objects lying on both red-sequences as cluster galaxies (Fig. [I]). 
This yields 258 early-type cluster galaxies down to K=23.4. We also include 9 galaxies which were not 
located within the red-sequence, but for which we measured a spectroscopic redshift equal to the redshift of 
the cluster (i.e. for 0.16 < z < 0.22, blue points on Fig. [I}. This cluster galaxy population is incorporated 
into the lens model as potentials described by Eq. [TJ with parameters scaled as a function of luminosity: 



1/2 



Tout 



1/2 



(TO 



1/4 



(2) 



The scaling relation for cr assumes that mass traces light, and its origin resides in the lFaber &: Ja ckson 
relation, that ha s been found to be reliable for describing early-type cluster galaxies (e.g. JWuvts et al.ll2004 ; 
Fritz et aDl2005h . Since the mass M scales as CTq r cut, we have: M oc L, assuming that the mass-to-light ratio 
is constant for all cluster members (note however that the mass-to-light ratio is not constant as a function 
of radius as we are finding that r cut > ineffective)- For a given L* luminosity, we will search for the values of 
<7q and r* ut , the only two free parameters, that gives the best fit, while r* orc is fixed at 0.15 kpc. Note that 
there are other possible scaling relations for r cu t, that we do not investigate in this work. 

The other parameters describing the small scale clumps associated with early-type galaxies are set as 
follows: the center of the dark matter halo is assumed to be the same as for the luminous component, 
and the ellipticity and position angle of t he mass is assumed to b e the same as the ones of the light. These 
assumptions are supported by the work by Koopmans et al. ( 20061 ) . who found that for a sample of early type 
lens galaxies: (i) isophotal and isodensity contours trace each other, and (ii) the position angle alignment 
between the stellar componen t and t he singular isothermal ellipsoid lens model are coincident. However, we 
caution that Koopmans et al. (2006) considered early- type field galaxies in their studies, and we extrapolate 
their results to early type galaxies residing in a very massive galaxy cluster. It is not clear whether or not 
this extrapolation is legitimate, and at this stage we cannot prove th is assu mption further. Galaxy-g alaxy 
lensi ng studies from the C OMBO-17 group (jKleinheinrich et al.l 120061 ). SDSS (jMandelbaum et al.l 120061 ) and 
RCS (jHoekstra et al]|2004l ) also suggest that the ellipticity of the light and mass are well correlated. But we 
caution once again that these galaxy-galaxy lensing studies were for field galaxies. 



Fig. 1. — Color magnitude diagrams and the selection of cluster member galaxies. The red-sequence selection 
is shown in the red boxes: all galaxies in this box are considered to be early-type cluster galaxies. The red 
points correspond to rejected galaxies, as we measured a redshift different from the cluster redshift. The 
blue points correspond to galaxies for which we measured a redshift equal to the cluster redshift, but that 
were not identified in both red sequences. 
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3.2. Source Plane Optimization 

We perform the optimization in the source plane, i.e. we optimize the fit by mapping the positions of 
the resolved multiple images back to the source plane and requiring them to have minimal scatter. 

The motivation for optimizing in the source plane is that it is much faster in terms of computing time 
than optimizing in the image plane. We are aware of the likely bias of magnification when optimizing in 
the source plane, but we verify that it becomes equivalent to fit in the source or the image plane when the 
number of multiply imaged systems increases (see H06). Beside, for images located close the the critical 
lines, the % 2 calculation in the image plane is sometime uncertain, making the source plane inversion more 
secure and stable. Furthermore, the total number of images that constitute a multiply imaged system may 
be unknown. Some of them can be demagnified and not detected, making image plane optimization delicate, 
whereas any observed image can be assigned a source position in the source plane. 

A x 2 estimator is constructed to quantify how well each trial lens model fits the observational data. 
Considering j multiply imaged systems and k critical curve constraints, we have: 

X 2 =EXpo S W+EXcrit(fc) (3) 

i k 

The first term, a sum over Xp OS (j), is constructed as follows: given a model, for each multiply imaged 
system j, we compute the source location uf = (xf,yf) for each observed image uj (1 < i < n), using the 
lens equation: uf = uj — Vy(u|), where ip is the lens potential. Then the barycenter u B — (x B ,y B ) is 
constructed from the N sources, and we define: 



^pos 



To estimate cr p0 s, we assumed that the typical uncertainty in measuring the position of any image is <7j = 0.2". 
This uncertainty is lensed back into the source plane, using the amplification A: erf = A cr| and we take as 

as (Tpos • 

The second term, a sum over x^ rit (fc), measures how well the locations of points on a critical line are 
reproduced by the model. This stringent constraint can be used in cases where we know with great confidence 
where the critical curve passes through. We define (x°^ s , y°^ s ) and (x™ od , y^t° d ) as the observed and model 
critical line positions respectively and construct XcritJ where Ax cr it and Ay CI i t are the uncertainties in the 
location of the critical lines (set to 0.3"): 

2 _ (^ s -xT d ) 2 +(y°t s -yT d ) 2 ^ 

Xcrit A 2 i A 2 vv 

This additional constraint in critical line position has been used in two cases (fc = 2): for the radial pair in 
system 5, we imposed the location of the radial critical line at z — 2.63 to bisect this radial pair; and for 
system 12, we impose the location of the tangential critical line at z = 1.83 to bisect the tangential pair. In 
practice, we have x^ os >> xlnt and X 2 ~ Xp OS - 

For such a complex cluster of galaxies as Abell 1689, the former parabolic optimization available in 
lenstool does not work well as it is too involving and time consuming. Furthermore, lenstool provided 
no estimate of the errors on the opti mized parameters, unless dedicated investigation of space parameter 
around the best model is performed ( Smith et al. 200a). In order to measure errors and to avoid local 
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X 2 minima, we implemented the Bayesian Monte Carlo Markov Chain (mcmc) package bayesys (Skilling, 
2004) in lenstool. According to a user defined model parameter space, the MCMC sampler draws randomly 
mock models and compute their % 2 . Progressively, it converges to the most likely parameter space and then 
it draws some realizations in this parameter space. We use these models to do statistics, in particular to 
compute error bars on the estimation of the parameters. With this new optimization method, lenstool 
returns samples of points but also the evidence associated to the model we are optimizing. The evidence 
characterizes our confidence in the model we optimize. It works as Occam's Razor (MacKay 1991) : if a 
simple model with few free parameters gives a good fit, a complex model giving the same fit will ha ve a 
lower evidence. We describe the Bayesian mcmc method in a dedicated publication (j Julio et al.ll2007l ). 

In the following, in order to quantify the goodness of a fit, we quote in Table [T] the RMS both in the 
source plane and in the image plane. Note that in this modeling we only consider the positions of the images 
or the critical curves as constraints, i.e. we do not include any constraints on the flux in this analysis. 
In principle, any optimization process should include flux constraints. However, the optimization of the 
positions is dominant, and the work by Zekser et al. ( 20061 ) has demonstrated that adding flux constraints 
did not have a strong impact on the properties of the final recovered mass model. However, some multiple 
images have shown unexplained surface brightness anomalies, e.g. system 12. That should deserve a peculiar 
study in term of caustic perturbation at very small scale. Also a supernova could affect the relative brightness 
distributions confusing the relative fluxes. 



Fig. 2. — Color image from F475W, F625W and F775W observations. North is up, east is left. The 
substructure associated with the north-east bright galaxy groups is clearly seen in this image. Size of the 
field of view is 160" x 160", corresponding to 485 kpc x 485 kpc. New spectroscopically confirmed multiply 
imaged systems proposed in this work are shown (i.e. new systems not proposed in previous works and for 
which at least one of the image has been observed spectroscopically). Note how system 31 shows an Einstein 
cross configuration. 



Fig. 3. — Multiple images considered in this work, both candidates (red squares) and spectroscopically 
confirmed (blue circles). The central galaxies as well as one in the north-east galaxy group have been 
subtracted for clarity. Size of the field of view is 160" x 160", corresponding to 485 kpc x 485 kpc. 
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Table 1. Multiply imaged systems considered in this work 



Id ra dec F775W 

-2phot a z modcl b z spcc c R.MS d (source plane, ) RMS d (image plane, ) 



1.1 13:11:26.447 -1:19:56.68 23.44 3.03^ 3.0 0.41 1.82 

1.2 13:11:26.287 -1:19:59.69 23.69 3.04±° J3 

1.3 13:11:29.771 -1:21:07.34 24.52 3.27+.°'5g 3.0 

1.4 13:11:33.055 -1:20:27.31 24.16 2.94±°'^ 

1.5 13:11:31.933 -1:20:05.69 24.62 3.35+.°;" 

1.6 13:11:29.851 -1:20:38.32 25.82 1.06 + Jj£ 



2.1 13:11:26.528 -1:19:55.08 23.37 2.62+JJ'^ 2.5 0.34 1.62 

2.2 13:11:32.961 -1:20:25.31 24.18 2.57Z° % 

2.3 13:11:31.973 -1:20:06.89 24.36 2S4Z° Ts 

2.4 13:11:29.811 -1:21:05.94 24.48 2.36+.°;J{ 2.5 

2.5 13:11:29.878 -1:20:39.32 25.63 1.59+°!® 



3.1 13:11:32.043 -1:20:27.41 26.65 5.48 3.25 ±1.47 0.32 0.55 

3.2 13:11:32.174 -1:20:33.31 26.85 5.45 t.°'ll 

3.3 13:11:31.680 -1:20:55.88 



4.1 13:11:32.174 -1:20:57.33 24.62 imZ" ) 27 1-1 0.42 1.29 

4.2 13:11:30.518 -1:21:11.94 23.91 1-32+1 J° 

4.3 13:11:30.759 -1:20:07.89 25.12 1.47^^ 

4.4 13:11:26.287 -1:20:35.11 24.67 1.33+.°'3i 

4.5 13:11:29.838 -1:20:29.11 



5.1 13:11:29.064 -1:20:48.33 24.42 3.29I055 2.6 0.72 2.11 

5.2 13:11:29.224 -1:20:43.92 24.92 3.16t°ll 

5.3 13:11:34.109 -1:20:20.90 25.26 2.15±°;^ 2.6 



6.1 13:11:30.759 -1:19:37.87 23.61 1.22+°2| 1.1 0.54 3.35 

6.2 13:11:33.335 -1:20:12.10 23.85 1.31+° f 1.1 

6.3 13:11:32.734 -1:19:54.28 23.02 0. 94+.°'25 

6.4 13:11:32.481 -1:19:58.69 24.00 1.09±°H 



7.1 13:11:25.446 -1:20:51.53 23.31 4.92Zg 78 4.8 0.50 1.07 

7.2 13:11:30.679 -1:20:13.70 24.19 5.20tog5; 4.8 

7.3 13:11:29.824 -1:20:24.71 28.23 0.77tn°\ 



8.1 13:11:32.300 -1:20:50.78 24.54 2.63 + °;j| 2.30±0.21 0.55 2.49 

8.2 13:11:31.393 -1:21:05.59 24.31 2.77+°J° 

8.3 13:11:31.500 -1:20:13.95 25.71 2.75+°gg 

8.4 13:11:25.520 -1:20:20.15 23.97 OJOZq'H 

8.5 13:11:30.325 -1:20:30.46 27.42 0.77 + „|| 



9.1 13:11:30.298 -1:19:48.33 25.92 4.97+°™ 2.69±0.27 0.36 1.60 

9.2 13:11:33.515 -1:20:50.38 27.53 1.06t°" 

9.3 13:11:28.723 -1:21:15.60 25.72 5.16+°|i 

9.4 13:11:26.267 -1:20:26.56 27.15 5.17+°|J 



10.1 13:11:33.976 -1:20:50.93 23.66 1.75+035 1.8 0.19 0.43 

10.2 13:11:28.049 -1:20:12.30 23.30 1.54+„'H 

10.3 13:11:29.317 -1:20:27.71 25.15 2.57+°;^ 



11.1 13:11:33.333 -1:21:06.74 24.03 2.91 + „ l\ 2.5 0.22 0.60 

11.2 13:11:29.064 -1:20:01.09 23.36 2.87+5J" 

11.3 13:11:29.491 -1:20:26.31 26.78 1.58°'n 3 r , 2 



12.2 13:11:27.355 -1:20:54.73 24.30 1.99 + £|| 1.8 0.04 1.54 

12.3 13:11:27.208 -1:20:51.53 23.97 1.99 + °;f| 



13.1 13:11:32.821 -1:19:24.11 24.10 1.02 + „„ 1.5±0.5 0.38 4.72 

13.2 13:11:32.994 -1:19:25.51 23.68 0.72+jJ;^} 

13.3 13:11:33.408 -1:19:31.11 24.12 l.lO+g^f 



14.1 13:11:29.024 -1:21:41.77 24.98 3.37 + JJ;|2 3.4 0.12 4.02 

14.2 13:11:29.454 -1:21:42.67 25.99 3.64+°" 



15.1 13:11:28.083 -1:20:14.75 25.69 1.99+°;jjj 1.8 0.15 0.38 

15.2 13:11:34.076 -1:20:51.38 25.79 2.00+°;|| 

15.3 13:11:29.244 -1:20:27.36 27.16 1-97+n^o 



16.1 13:11:27.989 -1:20:24.96 23.68 1.81+jj;" 1.76±0.16 0.47 2.78 



- 13 - 



Table 1 — Continued 



Id 


RA 




F775W 


Zphot" 


b 

Zmodcl 


^apcc 


RMS d (source piano, ") 


RMS (image plane, " 


16.2 
16.3 


13:11:28.923 
13:11:34.396 


-1:20:28.36 
-1:20:46.37 


25.06 
25.35 


2 26 + - 43 

1 on + 0.71 










17.1 
17.2 
17.3 


13:11:30.655 
13:11:30.388 
13:11:24.983 


-1:20:24.76 
-1:20:27.56 
-1:20:41.37 


24.19 
25.11 
24.31 


2 74+ - 49 
2 02+ ' 40 

o oc + 0.43 
Z ' ZO -0.43 




2.6 


0.55 


1.13 


18.1 
18.2 


13:11:28.243 
13:11:33.809 


-1:20:09.34 
-1:20:54.58 


25.00 


2 56 + - 47 
Z ' DD -0.47 




1.8 


0.07 


0.19 


18.3 


13:11:29.364 


-1:20:27.16 


26.78 


1 oS + °- 73 










19.1 
19.2 
19.3 
19.4 
19.5 


13:11:31.633 
13:11:25.253 
13:11:31.953 
13:11:32.033 
13:11:30.218 


-1:20:22.56 
-1:20:19.75 
-1:20:59.18 
-1:20:57.38 
-1:20:33.76 


24.83 
25.23 
25.88 
24.76 
27.94 


, 79 + 0.36 

1 - ' -0.36 

2 74+0- 49 

_ +0 .34 
-0.34 
, F i n 47 

2- 58+. °,; 47 
4 54+ 1 - 66 

*-°*-0.73 




2.6 


0.78 


1.75 


21.1 
21.2 
21.3 


13:11:31.019 
13:11:30.792 
13:11:25.253 


-1:20:45.77 
-1:20:44.57 
-1:20:10.95 


25.15 
26.70 
25.59 


I 7Q + 0-37 
i- 1 "-0.37 
1 >5q + 0-34 
^- J »_0.34 
i —,-. + 36 
A -' 8 -0.36 


1.91±0.19 




0.35 


1.76 


22.1 
22.2 
22.3 


13:11:29.694 
13:11:29.614 
13:11:32.404 


-1:20:08.54 
-1:20:23.56 
-1:21:16.00 


23.66 
25.49 
23.27 


1 QQ + °-39 
J-- 33 — 0.39 

1 99+ " 59 

0.39 

1 9fi+ - 39 
^ — 0.39 




1.7 


0.25 


0.69 


23.1 
23.2 
23.3 


13:11:29.524 
13:11:29.551 
13:11:32.661 


-1:20:09.74 
-1:20:22.76 
-1:21:15.20 


24.63 
25.99 
24.72 


2 n o+0.40 

1 99 +0 - 62 
9 no ^ 

z - uu _0.39 


1.78±0.20 




0.14 


0.44 


24.1 
24.2 
24.3 
24.4 


13:11:29.187 
13:11:32.057 
1 Q.I 1 .-in 9Q£ 

13:11:33.712 


-1:20:55.98 
-1:19:50.33 
-1 -1 Q- i i Q9 
-1:20:19.75 


25.24 
24.80 
24.33 
24.99 


r-, 1-^-1-0 4fi 

2 - 5o :o.46 

o 4 o + 0.45 
Z - 4J _0.45 
01 +0.50 
^•81_ 69 




2.6 


0.51 


4.52 


26.1 
26.2 
26.3 


13:11:25.146 
13:11:31.326 
13:11:30.245 


-1:20:32.36 
-1:20:25.16 
-1:20:32.36 


24.67 
25.39 
27.42 


i.08tS;3» 

1 04 4 " 027 
- L ' u *-0.27 
77 + 2 -53 
u ' ' ' -0.23 


2.33±0.29 




0.66 


2.15 


27.1 
27.2 

97 *i 


13:11:25.174 
13:11:31.366 

1 o.i i .on iQi 

±o. i ±.ou. ly i 


-1:20:33.11 
-1:20:24.56 

- L.ZVI.OZ. i O 


25.22 
26.19 

on CO 


i - Si -0.37 
„ .„ + n 48 

!-58to.t4 

A cr + 1-63 

4.55_ 73 


2.38±0.28 




0.78 


2.56 


28.1 
28.2 


13:11:28.298 
13:11:34.260 


-1:20:10.91 
-1:21:00.01 


27.20 
26.47 


1 17+„o„ 

- L '- L ' -0.29 

2.oo±°;« 


2.93±0.48 




0.57 


1.71 


29.1 


io.ii.nn o i n 

13:11:29.240 


-1:20:57.78 


25.97 


,7 + 0.57 
2 - 47 -0.46 






0.41 


2.65 


29.2 
29.3 
29.4 


13:11:30.028 
13:11:32.164 
13:11:33.618 


-1:19:33.92 
-1:19:52.53 
-1:20:20.55 


25.00 
24.80 
25.84 


„ .__i_o.58 

e.n+ - 46 
Z - OU -0.46 
q qc+0.57 
°- OO -0-57 




2.5 






30.1 
30.2 
30.3 


13:11:32.424 
13:11:33.188 
13:11:33.662 


-1:19:19.50 
-1:19:25.81 
-1:19:32.51 


25.91 
25.80 
25.73 


4 49+ D - 72 

^-^ y -0.72 

q oq + 0.56 

76 

3.30lS:l 6 6 




3.0 


0.61 


11.20 


31.1 


13:11:30.362 


-01:19:51.13 


23.77 






1.8 


0.34 


1.1 


31.2 


13:11:33.271 


-01:20:44.37 














31.3 


13:11:28.960 


-01:21:10.19 


25.37 












31.4 


13:11:26.504 


-01:20:21.75 


24.80 












32.1 


13:11:32.190 


-01:20:03.34 










0.18 


0.75 


32.2 


13:11:33.218 


-01:20:20.75 








3.0 






32.3 


13:11:29.587 


-01:21:02.49 














32.4 


13:11:29.801 


-01:20:43.17 














32.5 


13:11:26.600 


-01:19:57.60 








3.0 






33.1 


13:11:28.440 


-01:21:00.39 


23.66 






4.5 


0.51 


1.56 


33.2 


13:11:34.646 


-01:20:33.56 


26.62 













35.1 13:11:28.493 -01:20:59.38 

35.2 13:11:33.952 -01:20:33.36 

35.3 13:11:29.427 -01:20:33.76 



1.9 



0.48 



1.52 
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Tablc 1 — Continued 



Id 


RA 


dec F775W 


-Sphot z modol z spcc 


RMS d (source plane, ") 


RMS d (image plane, ") 




36.1 
36.2 


13:11:31.563 
13:11:31.683 


-01:19:45.73 
-01:19:47.13 


3.0 


0.08 


1.15 


40.1 
40.2 


13:11:30.269 
13:11:26.174 


-01:20:11.43 
-01:21:02.54 


2.5 


0.39 


0.99 





a Photomctric rcdshift, from B05 

b Redshift estimation infered from the mass model when all spectroscopically confirmed multiply imaged system have been included 
in the optimization. Errors bars quote 3cr confidence level 

c Spectroscopic redshift, see Paper I 

d The mean scatter is given for the whole system, not for each individual image composing a system 
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4. Multiply Imaged Systems 

4.1. Methodology 

In considering multiply imaged systems, we benefit from the work of B05 as a starting point and we 
have followed their notation hereafter. For details on the multiply imaged systems, see B05 or H06. In 
Appendix we give more details and illustrations on the new spectroscopically confirmed multiply imaged 
systems we consider. Since the redshift of the images is crucial in constructing an accurate mass model, 
we first considered systems for which we have measured an accurate spectroscopic redshift for at least one 
of the multiple images. For each system, we verified that the geometrical configuration was coherent and 
checked the colors and morphology of the different members before including a system in the optimization 
procedure. Our multiply imaged system catalog in general agrees well with the one from B05. However, we 
had to correct the identification for a few systems, and we propose new multiply imaged systems. All multiple 
images included in this work are shown on Fig. [3[ and their coordinates and redshifts are given in Table [T] 
If the availability of a great number of redshift is a strong advantage it is still challenging to detect reliably 
all the images of a given system, specially the central ones that falls behind the central bright galaxies. We 
preferred excluding rather than including the images for which we are not sure of the identification. 

The very first model was constructed using as constraints the images constituting system 1 (z — 3.05) 
and system 2 {z = 2.53). Even at this stage, a second large dark matter clump, associated with the galaxy 
group at north-east, was needed to accurately reproduce the location of the images of these two systems. 
Then we included the other spectroscopically confirmed multiply imaged systems. Before doing so, we use 
the current model to predict the location of the different images that constitute this new system. To do so, 
we use the image(s) for which we have measured a spectroscopic redshift as an input. This is useful to check 
if any new multiply imaged system agrees well or not with the current model. For example, this procedure 
indicated that system 12, as proposed by B05, was problematic. 

In order to better retrieve the locations of some images, we included 3 galaxies explicitly in the opti- 
mization procedure (this means that the parameters describing these galaxies were also considered as free 
parameters and optimized, as explained in section [33] instead of being scaled as a function of luminosity). 
Indeed, these galaxies are located very close to some multiple images: Galaxy 1 (a = 197.859, S — —1.332) 
is close to some images belonging to systems 1 and 2. Galaxy 2 (a = 197.886, S = —1.332) is a bright 
galaxy close to system 6. We also included the BCG to better retrieve the central images and the radial 
system 5. The inclusion of these extra free parameters was made possible since the number of constraints 
we have is large enough. The parameters of these individual galaxies arc presented in Table [3l keeping in 
mind that there may be a degeneracy between them and the cluster parametrization. In total, we have 33 
free parameters in the modeling. 

4.2. New Systems 

The new multiply imaged systems are shown on Fig. O We propose five new spectroscopically confirmed 
multiply imaged systems: system 32 at z — 3.0 (Fig.[l3|); system 33 at z — 4.5 (Fig. [14]); system 35 at z = 1.9 
(Fig. rrS]); system 36 at z = 3.0 (Fig. [TSJ) and system 40 at z = 2.5 (Fig. [T7J). We have split system 12 as 
proposed by B05 into two systems: systems 12 and 31 fFig. [T2"|) . We began by following the work of B05, but 
did not succeed in reproducing the 5 images they quote in their work with our current model. By studying 
the individual images carefully, we found that the morphology of the different images were inconsistent 
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with each other. We decided to restrict system 12 to two images, 12.2 and 12.3, forming a small highly 
amplified gravitational arc which was already bisected by the critical line at z = 1.83 (the redshift of this 
arc). Moreover, object 12.1 in B05 is also spectroscopically confirmed at z = 1.83, being part of an other 
system, we renamed it 31.1. This object has a very characteristic shape, constituted by two spots. The 
model predicted 4 counter images, one of them being strongly demagnified in the central region of the 
cluster. We were able to identify the three brighter counter images with the correct relative orientation of 
the two little spots (Fig. [T2|) . We then constructed system 31 from these 4 images, giving an Einstein cross 
like configuration (Fig.[2|). Note that H06 first proposed the splitting of system 12 from B05 into two distinct 
systems. 

We also provide a secure measurement for the redshift of system 10 (z = 1.83), that was used in previous 
studies by B05 and HOG but with a wrong spectroscopic redshift estimation (z = 1.37). We do not consider 
systems 20 and 25 as proposed by B05 in our analysis. For system 20, we expect a third image which is 
not seen, and the morphology and colors of these two images are somewhat different. Moreover, we took a 
spectrum for the two images and found that if their redshift is similar, the individual spectra are different 
(see Paper I for details). Concerning system 25, we used image 25.2 from B05 to define a new system that 
we call 33. 

H06 also exclud ed system 20 from their analysis and proposed the same identification as we do for 



system 12. Note that lZekser et al.l ([2006) also reported some problems with both systems 12 and 20, finding 



a considerably improved goodness of fit when removing them from their analysis. 



4.3. Candidate Systems 

At this stage of the modeling, we found that the mass of the second large scale dark matter clump was 
not well constrained, as suggested by the large error bars on the estimated parameters of the clump. We 
interpreted this as due to the lack of observational constraints in this region, as systems 6 and 30 are the 
only spectroscopically confirmed systems in this region included in the modeling. On the other hand, the 
position of the main dark matter clump was found to be very well constrained. 

We used the current model to predict the redshifts for all the other candidate multiply imaged systems. 
The results are given in Table [TJ where the errors bars quote the 3a confidence levels inferred from the 
Bayesian mcmc estimation. The estimates of the other candidate multiply imaged systems are rather well 
constrained by the mass model and agree well with the photometrically estimated redshifts. All the candidate 
systems are assigned a redshift as estimated by the mass model. These systems are then included in the 
optimization procedure as additional constraints. We thus take benefit from all the available information 
to further constrain the mass model. Note that even if the redshift of a multiply imaged system is not 
known exactly, it is worth using it in the optimization procedure since it will give us a net positive number of 
additional constraints. For a system with known redshift, and N multiple images, we get 2 (n — 1) constraints, 
but if the redshift is not known the number of constraints reduces to 2 (n — 1) — 1 which is positive as N > 1. 



5. Strong Galaxy-Galaxy Lensing Events at Large Cluster-Centric Radii 



Inside the critical region (R < 100" from the center), multiple images arise due to the cluster potential. 
Sometimes they appear close to individual galaxies that locally boost the lensing signal. 
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We also found some remarkable strong lensing features, located outside the critical region of the cluster 
(see Fig. 0]). Strong galaxy-galaxy lensing features have been observed around many isolated galaxies, in 
particular around early type galaxies in the field. They are still rare events, because only a very favorable 
alignment between source, lens and observer can produce these strong lensing events, due to the relatively 
small mass of a single galaxy. The typical density of such events is generally rare in galaxy field observations 



with about one strong galaxy-galaxy lensing event for about 1000 elliptical galaxies (jMiralda-Escude &: Lehar 



1992). In this work, we detected five (of which three are clear lensing events) strong galaxy-galaxy lensing 
events in a very small portion of the sky (~ 202 2 — ir x 100 2 sqr. arcsec. ~ 2.8 sqr. min., i.e. the size of 
the entire field minus the size of the critical region of the cluster). They would be very unlikely without 
the presence of the cluster potential, which adds a strong external lensing signal (7 ~ 0.22 and k ~ 0.36 at 
100" from the center) to the lensing signa l of the individ ual galaxies, making strong lensing possible around 



these individual galaxies (see the work by iKovnerl (|1987l ) on these so-called 'marginal lenses'). These simple 



observations suggest that this galaxy cluster is very massive. It would be interesting to compare the density 
of strong galaxy-galaxy lensing events that are found (or not found) for other clusters. This density could be 
used in principle to diagnose the cluster's surface mass density in this transition region between weak and 
strong lensing and bridge the two lensing regimes. 

These galaxy-galaxy lenses are of particular interest: with redshift measurements, we could put direct 
lensing constraints on the mass distribution of the lensing galaxies, stronger than the ones put on the BCG 
and Galaxies 1 and 2. This is because the influence of the cluster is not as strong as it is in the center of 
the cluster, and therefore there will be less degeneracy between the cluster and the galaxy parameterization. 
Moreover, with high re solution spectrosc o py, w e could also combine lensing analysis with velocity dispersion 
profiles as proposed by Koopmans et al. ( 200(jh . and test the hypothesis we introduced before for modeling 



the cluster members (section 3.1). We list some properties of these galaxies in Tabled 



Fig. 4. — Strong galaxy-galaxy lensing events detected outside the critical region of the cluster: full color 
images constructed from F475W, F625W and F775W observations (inverse colors), and images with the lens 
galaxy subtracted. Top row, from left to right: galaxies number 1, 2, 3. Bottom row, from left to right, 
galaxies number 4, 5. Size of each panel is 30kpcx30kpc. The first three systems are clear strong lensing 
events, whereas the lensing interpretation is not as clear for the two last events. Spectroscopic measurements 
are needed to alleviate the doubt. 



6. Mass Distribution from Strong Lensing 



The main advance we have made in the inner mass modeling compared to previous works is to have 
included 34 multiply imaged systems as constraints, of which 24 do have spectroscopic confirmation. The 
positions of the observed images are reproduced accurately by our model for the majority of the images. 
The mean scatter in the source plane is equ al to 0.45". This is much better than the mean scatter in the 
source plane reported bv lZekser et alJ (|2006h of 1.19". In the image plane, the RMS is equal to 2.87" This is 
comparable to previous studies performing the optimization in the image plane: B05 reported a RMS equal 
to 3.2", and H06 2.7". As a complementary confirmation of the quality of the lensing model, we checked that 
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all spectroscopically confirmed background sources for which we did not identify multiply imaged systems 
were predicted to be singly imaged by our model. 

We begin this section by presenting the best fit PIEMD parameters we get as an output of the optimization 
procedurcQ In Appendix, we provide a comparison of the model presented here with other strong lensing 
works. We also discuss the interpretation of our model and the degeneracies we encountered. 



6.1. Parameters of the Dark Matter Halos 

The parameters found for the two large scale dark matter (dm) clumps are given in Table [3J The 
position of the main dark matter clump is found to be offset from the BCG by ~ 25kpc (Fig. O. This 
position is exactly where the lines connecting the two main radial arcs (systems 5 and 22) bisects. As 
we could have expected, this position is also where the line connecting images 31.1 and 31.3 and the line 
connecting 31.2 and 31.4 bisect (remember system 31 displays an Einstein cross configuration). The peak 
of the total mass distribution (corresponding to the superposition of the main dark matter clump and the 
BCG) however coincide with the BCG (green point on Fig. [5J. Thus the peak of the total m ass distribution is 



coinc ident with the x-ray center which is also found to coincide with the center of the BCG ([Andersson et al 



2006) 



The central mass distribution dominates the entire cluster, as suggested by the centrally peaked X- 
ray emission. However the second clump is also a substantial contributor and is really required by the 
observational constraints. This reinforces the idea that violent processes are ongoing in Abell 1689. 

We optimized over the parameters for three early-type cluster galaxies individually, rather than scaling 
their parameters with luminosity. This was necessary to explain observed multiple image configurations in 
their vicinity, as they were found to dominate the lensing signal locally. Thus, additional constraints were 



2 A parameter file containing all the following information, and which can be used with the publicly available LENSTOOL 
software, is available at http://www.dark-cosmology.dk/archive/A1689/ This file can be useful for making model based pre- 
dictions, e.g. counter-images of a multiple image candidate, amplification and mass map and location of the critical lines at a 
given redshift, and it will be updated. We also provide the mass map generated from the best fit model as a FITS file for direct 
use. 



number 


RA 


DEC 


dBCG 


L6I1S ^spec 


Lens magnitude 


R, 


1 


197.872 


-1.309 


114 " 


0.1758 


18.25±0.13 


2.0±0.1 


2 


197.897 


-1.345 


90 " 


0.1844 


18.58±0.16 


2.2±0.2 


3 


197.884 


-1.369 


110 " 


0.1855 


18.11±0.01 


4.1±0.1 


4 


197.894 


-1.323 


101 " 




18.69±0.15 


2.6±0.2 


5 


197.868 


-1.312 


105 " 


0.1990 


18.64±0.09 


3.0±0.2 



Table 2: Strong galaxy-galaxy lensing events detected outside the critical region of the cluster: coordinates 
(J2000), distance to the BCG; spectroscopic redshift measurements when available; F775W AB magnitudes 
obtained from the surface brightness profile fitting by H06; circularized physical half light radius in units of 
kpc (H06). 



- 19 - 



derived for these individual galaxies. We present these galaxies, and their best fit parameters in Table [3] 
One should however bear in mind, that the parameter constraints for these individual galaxies are likely 
somewhat degenerate with the parameters of the large scale dark matter clumps. This is particularly true 
for the clump we associated with the BCG and with Galaxy 2: these values are too high to be representative 
of galaxies. 

We included all the other identified cluster members, most of which were elliptical galaxies, in the 
optimization. We scaled the parameters as a function of luminosity using the scaling relation discussed in 
section 13.11 We considered a cluster galaxy with luminosity L* corresponding to a F775W magnitude of 
17.54. At the end of the optimization procedure, we get the best parameters describing this L* galaxy. We 
find Oq = 129.7 ± 3.3 kms -1 and r* ut = 53.5 ± 5.0 kpc (1 a). Note that there can still be some degeneracy 
between the mass we put in the large scale dark matter clumps and the mass associated with the galaxy 
cluster members. Keeping that in mind, we find galaxy halos to have a rather small spatial extent (compared 
to field galaxy halos). Thus we provide evidence for truncation of galaxy dark matter halos in high density 
environments. Even if we are aware of the possible degeneracy mentioned above, we are confident in these 
values since they are in good qualitative agreement with independent weak galaxy-galaxy lensing results and 
with the expectations from numerical simulations. 



6.2. The Cluster Mass Profile 



In Fig. O we show the R band image of Abell 1689 along with the contours generated from the projected 
mass map inferred from the best-fit model (red contours). These contours show where the projected mass is 
equal to 1.6, 2.4, 4.0 x 10 10 M^arc sec" 2 . This mass map is found to be in very good agreement with the 
presented by Zekser et al. ( 20061 ). By integrating this two dimensional mass map, we get the total mass 



one 



profile shown in Fig. [5] The horizontal line at the top of Fig. [5] shows at what distance, from the center 
of the cluster, the multiple images used to constrain the model are located, hence where the mass estimate 
is robust. In Fig. [6j we compare the mass profile derived in this work with the one from B05 and H06, 
and we find a very good agreement. As a result the inferred concentration parameter we find by fitting the 
inner mass profile with a nfw profile (C200 = 6.0 ± 0.6 at 3a confidence level) is in good agreement with 
previous estimations. This is encouraging since different studies have used di fferent techniques and d i fferent 



sets of multiply imaged systems. Also shown is the mass profile derived by lAndersson &: Madeiskil (|2004f ) 



from X-ray data. The observed discrepancy in these independent estimates will be discussed in section [8l 

To quantify the accuracy of the mass map, we considered 1 000 realizations coming out from the MCMC 
optimization procedure, and constructed a mass map for each one. The maps where then stacked and 
the mean and the standard deviation (dispersion) were computed at each point. Then this dispersion was 
converted in a percentage on the accuracy of the mass measurement. The accuracy of the mass reconstruction 
was found to be smaller than 15% over the whole field. A similar procedure was used to estimate the error 
bars of the mass profile presented in Fig. [6l 

Fig. 5. — R band image of Abell 1689. Size of the field of view is 160" x 160", corresponding to 485 kpc 
x 485 kpc. The red contours show where the projected mass density equals 1.6, 2.4, 4.0 10 10 M arcsec" 2 . 
The precision on the mass measurement, as inferred from the MCMC optimization, is found to be smaller 
than 15% over the whole field. The green blob shows where the peak of the mass map is found. The blue 
circle shows the location of the main dark matter clump (3 a limits) . 
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Fig. 6. — Total projected mass as a function of the aperture radius. The shaded area corresponds to the mass 
profile derived from this work, within the la errors as inferred from the mcmc optimization. The horizontal 
line at the top of the plot shows where multiple images are found, hence where the mass estimate is more 
robust. For comparison we plot the mass profiles from H06 (solid line), B05 (dashed line) and Andersson & 
Madejski (2004) from an x ray analysis (dotted line). For clarity, some error bars are excluded. 

6.3. Critical Lines 

The critical curves are well defined and illustrated in Fig. [7| for a source redshift equal to 3. The outer 
line corresponds to the tangential critical curve, whereas the inner line corresponds to the radial critical 
curve. 

Caustic curves correspond to the unlensed location of the critical curves. For a given system, the multiple 
image configuration will depend on the location of the source with respect to the caustics. Fig. [5] shows the 
radial and tangential caustics for a source redshift equal to 1.83 (redshift of system 12, 31 and 10) and the 
reconstructed location of the different sources (only the sources corresponding to spectroscopically confirmed 
multiply imaged systems are shown). The tangential caustic curve crosses source number 12, resulting in a 
highly amplified tangential arc in the image plane. 

It is interesting to look at the distribution of the strongly lensed sources. From Fig. [SJ we clearly see 
a group of 3 galaxies at z = 1.83 (sources 10, 12 and 31) that have been strongly lensed by the cluster. It 
has a maximal spatial extension of 6.5", which translates into 55 kpc. Moreover, the redshift differences 
between these three galaxies translate into velocity differences smaller than 300 kms" 1 . This means that 
these 3 galaxies are likely to be bound and interacting. Moreover, we have other candidate multiply imaged 
systems (systems 15 and 18) whose estimated redshift is compatible with z = 1.83. Thus it is possible that 
we are observing a galaxy group strongly lensed by Abell 1689. Interestingly, this galaxy group would stand 
in the so called 'redshift desert'. In addition we detect two pairs of likely interacting galaxies at z = 2.63 
and z = 3.05. The properties of these galaxies are described in Paper I. 



6.4. Mass Associated with the Galaxies 

We want to compare the mass associated with the individual galaxies (M ga i ax ), to the total mass, (M to t), 
as a function of radius (Fig. [9]). Inside the Einstein radius (Re ~ 136 kpc), we find M ga i ax = 13% M tot , 
whereas inside a radius of 364 kpc (roughly the ACS field of view), we find that the contribution of the galaxies 
to the total mass is 7%. As M ga i ax is proportional to the galaxy luminosity, these means that we have a 
smaller mass-to-light ratio in the central part of the cluster compared to the outer region, clearly showing th e 
effect of mass segregation in cluster center. Similar trend was also observed in C10024 (|Kneib et alJl2003h . 



We caution the reader about the possible degeneracy between the smooth dark matter component and the 
galaxy component. Moreover, we did not investigate other models to describe the galaxies. However, we are 
confident in the PIEMD results presented her e for the galaxy com ponent since they agree with independent 



studies of galaxy-galaxy lensing in clusters. IZekser et al.l ([20061 ) also compared the relative contributions 



Fig. 7. — Critical lines for a source redshift equal to 3. The outer line is the tangential critical line, the inner 
line is the radial one. Size of the field of view is 160" x 160", corresponding to 485 kpc x 485 kpc. 
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Fig. 8. — Radial and tangential caustics at z = 1.83 and the location of the reconstructed sources. The 
radial caustic is the one with a circular shape. 

from the galaxies and the halo, and found the galaxy component to account for 15% of the total mass within 
the Einstein radius. 

Fig. 9. — Contribution of the galaxy component to the total mass as a function of radius. The vertical 
dotted line shows the location of the Einstein radius Re (for a source redshift equal to 3). 



7. Mass Distribution from Weak Lensing 

As an additional test for our strong lensing model, we check the agreement with the weak lensing analysis 
at large radius. Previous works have found a discrepancy between the NFW concentration parameters derived 
from weak and strong lensing separately. 



7.1. Shear Profile Measured from CFH12k Data 



We directly measure the shear profile of Abell 1689 in wide field CFh12k data. A shear profile is 
constructed from the shapes of the lensed background population, so any contamination of the background 
catalog by foreground members will dilute the shear signal. Any weak lensing study will need to minimize 
this contamination, although it is not always easy and often impossible to avoid for single band observations. 
However, this dilution of the weak lensing signal by cluster members can also be turned to our advantage 
and use d to derive the properti es of the cluster population, in particular the radial light profile (see recent 
work bv iMedezinski et al.l l2006). 



Shear profiles for Abell 1689 have previously been measured usin g the same data set. Bardeau et al 



([20051 ) used a magnitude cut-off to select background sources, while iBardeau et al.l (|2007l ) used a color 
selection (defined by r-i > 0.7, corresponding to the color of the red-sequence members). The latter provided 
a less contaminated background catalog than the former, but the shear profile was still found to be shallow 
in the center, where the cluster contamination is more prominent. If the virial mass derived from these 
different studies roughly agrees, the estimation of the concentration parameter differ significantly. Below we 
will investigate in detail this contamination problem through several rejection criterion. 



Construction of the Background Catalog: We performed a Bayesian photometric study to estimate a 
redshift probability distribution (Pwss) for each object, based on the B, R, I photometry. The details of the 
method are given in iLimousin et alJ (|2007l ). where we showed that this photometric redshift estimate was 
reliable for galaxies with z > 0.5. We then apply a criterion to determine whether a galaxy is a foreground 
or a background galaxy, setting the limit at z — 0.3. In particular, we first compute the quantity : 



1 f +oc 

N J z 



(6) 
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where N = / +o ° Pba, yC s{z')dz' is a normalization factor. This quantity provides information on the fraction 
of the probability distribution which is beyond z. In practice, the upper bound of the integral is set to 5. 

We search for the threshold value of z, which best discriminates the cluster members, while still leaving 
enough background galaxies for a weak lensing study. To do s o, we use a su b-sample of 239 elliptical cluster 
galaxies for which we have spectroscopic redshifts (data from ICzoskd 120041 ) . These objects are located out 
to 340" from the center of the cluster, i.e. where the contamination is most likely to be strong. We found 
that if we impose Xo.4 > 60, we actually get rid of most of the contamination by the cluster members (with 
only 4% contamination remaining), while at the same time keeping enough objects to pursue a weak lensing 
study, with a background source density of 12 galaxies/arcmin 2 . 

The spectroscopic subsample of red cluster members used to calibrate our rejection criterion does not 
allow to test the removal of the faint blue cluster members which are also a source of contamination at large 
radius from the cluster center. Considering the R-i color of the objects which fulfill the criterion xo.4 > 
60, the majority (~ 60%) are red objects, with r-i > 0.7. This means that our criterion does select blue 
objects that could be also cluster members. To test this residual contamination, we compared the shear 
profile presented in Fig. [TO] with the one inferred from a catalog whose members fulfill both \oa > 60 and 
R-i > 0.7. We find both shear profiles to be comparable, in particular in the inner part. We are therefore 
confident that our selection criterion does not include too many blue cluster members. 

The final catalog used to construct the shear profil e contains objects for whi ch Yn.4 > 60. F ig L [TOl gives 
our shear profile, along with the shear profiles derived by lBardeau et al.l (|2005l ) and lBardeau et al.1 (|2007l ) . All 
shear profiles agree well beyond R ~ 400", where the contamination by cluster members becomes negligible. 
Closer to the center of the cluster, the differences between the different shear profiles increase, suggesting 
that previous analysis do not minimize the contamination. 



Discussion: Although not perfect, the rejection criterion we propose in this work seems to exclude most 
of the cluster members, as it gives an expected non-flat shear profile, in agreement with the predictions from 
the strong lensing model (see next subsection). Another advantage of this rejection method is that it allows 
keeping the distant lensed objects whose colors are as red as the cluster red sequence, whereas a color cut-off 
will reject these objects. In conclusion, we are able to keep more background objects than when applying a 
color-cut off (10 280 instead of 6 280), resulting in a shear profile with higher signal to noise ratio. 



The analysis by iLeonard et al.l (|2007l ) provides an estimation of the tangential reduced shear of about 
0.2±0.03 at a distance of 100", in good agreement with our measure ments (Fig. [TOl) . This is encouraging for 
the rejection criteria we propose, especially because the analysis bv lLeonard et al.l ( 2007 ) relies on the ACS 
data and is likely not to suffer from contamination problems as can be the case with wide field ground based 
studies. 

Nevertheless, this improved rejection criterion described here keeps some limitations. From the spec- 
troscopic sub-sample, we estimate a 4% remaining contamination, but this may be underestimated since the 
spectroscopic sub-sample is not fully representative of the whole cluster population. Moreov er, our Bayesian 



photo metric redshift analysis is based on only three colors informations. We have shown (jLimousin et al 



20071 ) that the redshift estimation was correct, but ultimately this background galaxy population should be 
constructed from a more robust photometric determination, using many filters. 
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7.2. Shear Profile Predicted by the Strong Lensing Model 

The strong lensing analysis shows that the cluster has two large scale dark matter clumps which dominate 
the mass distribut ion of the cluster a s a whole, and no significant clump is found outside the region probed 



by the ACS data (jBardeau et al.ll2005l ). In addition, the x-ray emission is peaked near the center of the mass, 
making us confident that we have constrained most of the mass of the cluster. We should therefore be able 
to provide reliable predictions for the properties on large scales, although they should be strictly considered 
as lower limits. 

We consider a background population distributed in a plane at z ~ 0.9, which is the mean estimate 
of the Bayesian photometric redshifts for the background catalog we constructed. Background galaxies are 
randomly distributed on the plane, with a density of 40 galaxies per square arcminute and shapes are assigned 
by drawing the ellipticity from a Gaussian distribution similar to the observed CFHT ellipticity distribution. 
We only consider the intrinsic noise corresponding to the width of the ellipticity distribution. In particular, 
we assume that the PSF can be perfectly estimated and subtracted and do not introduce any additional 
observational noise. We calculate the deformations produced by our best-fit mass model for this cluster on 
this background galaxy population. The resulting catalog of lensed images is then used to construct a shear 
profile out to large radius. This shear profile is plotted as a solid line in Fig. [TO] and is referred to as the 
predicted shear profile. 

We are aware of the limitations of the predictions we are making here: the strong lensing constraints 
span up to 85" from the cluster center, thus the predictions on larger scales will be model dependent. The 
farther away from the center of the cluster, the less reliable the predictions from the strong lensing model, 
and the only meaningful comparison between the weak and strong lensing regime can be made in the region 
of overlap between the data sets. Strictly speaking, only the three first data points of the measured reduced 
shear profile on Fig. \W\ can be compared to the predictions of the strong lensing model (i.e. up to 150" 
from the cluster center). Note that this inner region is precisely where the disagreement between shear 
measurements is the largest. The method exposed in this subsection can be used to predict the level of the 
shear signal which is expected and compare it to the measured weak lensing signal. This can be useful to 
test for contamination in weak lensing catalogs. Fig. [TOl shows that the shear profile predicted by the strong 
lensing model is found to be in agreement with the shear profile measured using weak lensing. 



The shear profile measured by iBroadhurst et al.l (|2005bl ) using SUBARU data, is shown in Fig. [TTJ and 



compared to the one presented in this work, based on CFHT data. We also show the curve of the best NFW 
fit to our weak lensing measurements. This fit is found to be good. Comparing SUBARU and CFHT reduced 
shear profiles, we find marginal agreement. In particular, in the central 150", the SUBARU shear profile is 
higher than the CFHT one, reaching values as high as 0.7 in the very center. As a consequence, when fitting 
an NFW profile to the S UBARU shear profi l e only , high values for the concentration parameter are derived: 



H06 finds C200 ~ 30, and Medczinski et al. ( 20061 ) finds C200 ~ 22 when using a slightly different background 



galaxy catalog than IBroadhurst et al.l (|2005bl ). This will be discussed in the next subsection. 



Fig. 10. — Comparison between shear profiles constructed using different rejection criteria to select back- 
ground lensed sources: black squares - corresponds to the Bayesian photometric redshift based selection; red 
circles - a color selection, and blue crosses - a magnitude cut-off is used (some of the errors bars omitted for 
clarity). The number of objects corresponding to each rejection criteria is given in brackets. We find that 
all shear profiles agree well beyond R ~ 400". The solid line corresponds to the strong lensing prediction, 
and the green filled triangle comes from the study by Leonard et al., 2007. 



Fig. 11. — Comparison of our measured shear profile (black squares) with the one measured from SUBARU 
data presented in Broadhurst et al. 2005b (red circles). The dashed line corresponds to the best NFW fit to 
our shear profile, and the green filled triangle comes from the analysis by Leonard et al., 2007. 



7.3. The NFW Concentration Parameter 



The different studies performed on Abell 1689 have given discrepant values for the concentration pa- 
rameter. In our study, we find the strong lensing to agree well with the weak lensing analysis and gives 
C200 ~ 7. 



Fitting an NFW profile to the reduced shear profile of Abell 1689 (Fig. [TTj) , we find: r2oo = 2.16 ± 0.10 
Mpc and C2on = 7.6 ±1.6. Considering the background catalog constructed with a magnitude cut-off, 
Bardeau et all |2005h found: r 2 on = 1-9 Mpc, C onn = 3.51q'3, and considering the background catalog 
constructed with a color selection. iBardeau et al.1 (|2007l ) found: r 2 ao — 2.2 Mpc, c 2 oo = 5.2+° '|. Comparing 



these values, we see that T2qo (and hence M200) is a reliable estimator that does not depend strongly on 
the contamination of the background catalog by cluster members. This shows that the mass estimate from 
weak lensing is a robust quantity, even if they are based on contaminated background catalogs (as is often 
the case, in particular for studies with data in one band only). On the other hand, the values inferred for 
the concentration parameter strongly depend on the inner part of the shear profile which is more affected by 
the contamination. As less contamination in the background catalog means a steeper shear profile, it leads 
to a higher estimate of the concentration parameter. 

It is interesting to compare the concentr ation parameters found in other studies since the values for 
r 2 oo roughly agree with each other (Tabled]). IClowe fc Schneider! (|200ll ). considering a background catalog 
constructed with a magnitude cut-off, found: C2 on ~ 6, and when us ing a more re fined maximum likelihood 
method to fit the data, C200 = 4.7 was found king et al.ll2002bh . Iciowe! (|2003h . using color information 
from CFHT (the same as the one used in this p aper) to better s e lect ba ckground sources, and using a 



maximum likelihood method, found C200 ~ 7.9. iBroadhurst et al.l (|2005bh found a higher concentration, 
c v ir ~ 13 (corresponding to 0200= 10.8) when fitting simultaneously both strong and weak lensing data 
points. However, H06, when fittin g their strong lens i ng bas ed mass profile with the weak lensing SUBARU 
data, found C200 = 7.6jlQ*J, when IBroadhurst et al. ( 2005b[ ) was finding c 2 oo ~ 10.8. This is surprising, 
because the inner mass profiles inferred by B05 and H06 from strong lensing agree well, and both studies 
use the same weak lensing data from SUBARU. The only difference between these two fitting is that B05 fit 
the NFW profile to the surface mass density map obtained from the shear values, whereas H06 fit the shear 
directly. But this is not likely to explain such a large difference, since both quantities are two different ways 
of characterizing the same mass distribution. The only explanation may come from the relative weights 
that have been assigned to the weak lensing and strong lensing data points when constructing an overall \ 2 
since both regimes are found in disagreement by B05 and H06. Assigning different weights to the different 
data points will lead the x 2 to be more representative of the strong lensing or the weak lensing data points. 
Another contribution may come from the different errors estimations in the strong lensing mass profile found 
by B05 and H06. This demonstrates that if the weak lensing and strong lensing regime do not match each 
other, it is not possible to give a reliable estimation of the concentration parameter. In that sense, the value 
of 7.6 found by H06 is misleading because this provides a bad fit to the SUBARU weak lensing data (see 
Fig. 21 of H06 and compare with our Fig. [TT|) . over predicting the outer radial profile. As a consequence, 
they find r 2 oo = 2.55 Mpc, which is larger than the value found in this work, and which leads to a M 2 oo 
estimate which is twice the one presented here. In this work the fact that the strong lensing profile and 
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Table 3. piemd parameters inferred for dark matter clumps considered in the optimization procedure. 
Coordinates are given in arcseconds with respect to the BCG. The ellipticity e is the one of the mass 
distribution, expressed as a 2 — b 2 /a 2 + b 2 . Error bars correspond to la confidence level as inferred from the 
mcmc optimization. Values into brackets are not optimized. 



Clump 


R A 


DEC 


e 





fcore (kpc) 


r cut (kpc) 


cto (kms ) 


Clump 1 


0.5 ± 0.2 


-8.5 ± 0.4 


0.21 ±0.01 


90.4±1.0 


99.4±3.9 


[1500] 


1212.7±13.0 


Clump 2 


-70.7±1.4 


49.1±3.4 


0.7±0.03 


78±2.5 


66.6±6.0 


[500] 


542.8±3.6 


BCG 


-1.0±0.2 


0.1±0.4 


0.47±0.04 


66.0±6.3 


5.2±1.0 


128.5±37.0 


370.5±10.3 


Galaxy 1 


[49.0] 


[31.5] 


0.70±0.13 


113.5±9.0 


25.1±2.8 


158.3±16.5 


220.5±11.4 


Galaxy 2 


-45.6±0.5 


31.6±0.8 


0.77±0.04 


45.5±2.3 


17.2±2.2 


179.5±7.5 


357±20.2 


L* elliptical galaxy 










[0.15] 


53.5±5.0 


129.7±3.3 



C200 


Method 


Reference 


Remark 


G 
4.7 

7.9 

o c+0.5 
°-°-0.3 

5.2 ±0.3 
22.1^-9 


WL 
WL 
WL 
WL 
WL 
WL 


Clowe & Schneider (2001) 
Kinsr et al. ( 2002b) 
Clowe (2003) 
Bardeau et al. (2005) 
Bardeau et al. (2007) 
Medezinski et al. (2006) 


magnitude cut-off 
same data as Clowe & Schneider (2001) 
color selection 
magnitude cut-off (~ 25 000 objects) 
color selection (~ 6 300 objects) 
based on SUBARU data 


r 7 +0.34 
°- ' -0.5 

6 ±0.5 


SL 
SL 
SL 


B05 

Zekser et al. (2006) 
H06 


non-parametric method 
nfw (15%) + Shapelets (85%) 
parametric method 


7 fi +0 - 3 

' -°-0.5 


WL + SL 
WL + SL 


Broadhurst et al. (2005b) 
H06 


ACS + SUBARU data (fitting surface mass density) 
ACS + SUBARU data (fitting shear profile) 


7 7+ 1 - 7 

'■'-2.6 


X-ray 


Andersson & Madeiski (2004) 


hydrostatic equilibrium assumption 


7.6 ± 1.6 (lcr) 
6.0 ±0.6 (3a) 


WL 
SL 


This work 
This work 


bpz selection (~ 10 300 objects) 



Table 4: Concentration parameters found in different studies. 
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the weak lensing profile are fully compatible offer the best opportunity to derive a good estimation of the 
concentration parameter. 

The concentration parameters derived from our weak lensing analysis (C200 = 7.6 ± 1.6) agrees with 
the one derived from our strong lensing analysis (C200 — 6.0 ± 0.6). Thus we reconcile the concentration 
parameters derived from the two lensing regimes. This values also agree with the one inferred from X-ray 
data, where lAndersson fe Madeiskil (|2004 ) found c 200 



7 7+ 1 - 7 

' • ' -2.6- 



For the time being, the weak lensing studies based on wide field multi-color observations are giving very 
differ ent results for the c oncent ration parameter. Fitting an nfw profile to t he SUBARU weak lensing data only 
from lBroadhurst et alj |2005bl ) gives c 200 ~ 30 (H06), or c 20 o = 22.1^ (jMedezinski et alj|2006h . whereas 



fitting an nfw profile to the CFHT data presented here gives c 2 oo = 7.6 ± 1.6. The main difference between 
SUBARU and CFHT reduced shear profiles is in the central part (Fig. [TT|): the SUBARU reduced shear profile 
is higher, leading to higher concentration parameter than the CFHT reduced shear profile. Discriminating 
between high (c 2 oo > 20) and low (c 2 oo < 10) values of the concentration parameter is of interest since a 
high value of the concentration parameter, if confirmed independently, will have important consequences for 
cosmological models. Thus the weak lensing issue on Abell 1689 is not settled yet, and an independent study 
is needed to disc riminate betwe e n SUB ARU and CFHT shear profiles. Recently, an independent study provided 
one data point: iLeonard et al.l (|2007l ) measured the reduced shear to be equal to 0.2±0.03 at a distance of 
100" from the cluster center (Fig. [TT]) . in agreement with the CFHT data. This estimation is likely to be not 
contaminated by cluster members since the analysis relies on a small and deep ACS frame where the authors 
can control much better (from visual inspection and from the extensive spectroscopy available for central 
cluster members) the eventual contamination by foreground members compared to wide fiel d weak lensing 



studie s dealing with large catalogs of galaxies. Moreover, the deepness of the ACS data allow ILeonard et al 
([20071 ) to reach a density of background sources equal to ~ 200 sources/arcmin 2 . This unique data point can 



help to resolve the discrepancy between SUBARU and CFHT shear profiles reduced shears values at a distance 
of 100" from the center of the cluster. However, this data point is not enough to fully settle the point and 
another independent study providing a measurement of the reduced shear profile on a larger range of radius 
is needed. 



8. Discussion and Conclusions 



Fig. [6] shows the excellent agreement found by different strong lensing studies for the inner mass profile, 
though they use different methods. One of the main conclusion of this work is that a good agreement for 
the concentration parameter is found by the different methods used (strong and weak lensing and x-ray). 
We can see from Table [4] that the mean value of c 2 oo is around 7. Note that the X-ray measurement appears 
to be quite efficient to determine the concentration parameter. This is not surprising since so far it seems 
that this is the best way to determine the universal profile of dark matter halos, even providing one of the 
best proof of this universal profile. 

This value (c 2 oo ~ 7) is slightly higher than the ave rage value predicted b y N-body simulations of cluster 
formation and evolution in the Acdm cosmogony by iBullock et alJ (|200lh that predicts a concentration 
parameter c 2 oo ~ 5.5. It is possible to reconcile observed high values of the concentration parameter with 
the predicted one, e.g. invoking tri axiality of the dark matter halo ( Clowe et al. 2004 : Qguri et al. 20051 : 
GavazzilEool ICorless fc Kmell2006h . neighboring massive structures |King fc Corlessl |2007| ). or barvonic 
physics (which leads to more concentrated dark matter halo profiles as gas cools in their inner region, see 
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e.q. Gnedin et al . 2004). Indeed, none of the earlier simulations bv lBullock et al.l (|200lh contain a statistically 
significant sample of massive clusters that could likely produce the plethora of stro ng lensing phenomena 
displayed by clusters like Abell 1689. The Millenium simulation (jSpringel et al.ll2006l ). with a large box size 
(500 h" 1 Mpc) on a side would offer an appropriate comparison sample. What is needed at the present time 
is a distribution of best-fit concentration parameters for observed lensing clusters. Given the large diversity 
in observed massive, lensing clusters an adequate comparison can be made with simulations only with a large 
sample (20 clusters or so). 

This apparent disagreement between observed a nd predicte d concentration parameters may not be 
interpreted as a sign of failure of the CDM paradigm. Engl (|2000h has shown that the quality of the NFW 
fitting depends on whether the halo is in equilibrium, and that substructures degrade the fitting quality, 
because the NFW profile was found for equilibrium halos. Keeping in mind that Abell 1689 is composed by 
two massive large s cale dark mat t er hal os, this means that we have to discuss the nfw fitting parameters with 
care. The work bv lZekser et al.l (|2006[ ) also points out that t he NFW profile may not be a good description 
of the central part of Abell 1689. Recently, Saha et al. ( 20061 ) deprojected the mass map of Abell 1689 and 
obtained inner profiles consistent with p ~ r _1 , supporting the NFW predictions. However, they deprojected 
the mass map assuming spherical symmetry, which is a strong prior that may bias the measurement. 

It is worth noting a remaining discrepancy between the mass estimates from lensing and X-ray measure- 
ments. In general, x-ray estimates agree well with gravitational lensing estimates for clusters with a high 
concentration of central x-ray emis sion (i.e. rel axed cooling flow clusters) , but are seemingly in disagreement 
for less centrally peaked clusters ( Allen 1998h . Abell 1689, as discussed before, seems not to be fully re 



laxed and does not present any cool core, though the x-ray emission is clearly unimodal. Rasia et al. ( 2006h 
studied a set of five galaxy clusters, resolved at high resolution, in a hydrodynamic simulation, examining 
the systematics affecting the x-ray mass estimates. They showed that for a cluster undergoing a merger, 
the assumption of hydro-dynamical equilibrium led to the mass being underestimated by 30%, and that a (3 
model gave even more discrepant results. However, a 30% correction would still not be enough to reconcile 
the X-ray mass estimate with that from lensing for Abell 1689, but the discrepancy then becomes smaller, 
of the order of 40%. This remaining discrepancy may be understood if Abell 1689 is undergoing a merger 
along the line of sight, which remains a possible sce nario at this stage as suggest ed by the broad redshift 
distribution of the cluster members and the SZ study. Andersson fc Madejskil (2004) argue that if Abell 1689 
was undergoing a merger along the line of sight, the x-ray analysis would underestimate the total mass by a 
factor of 2. The simple model (two perfectly spherical clusters aligned exactly along the line of sight) used to 
derive the estimate is implausible, with a more realistic model leading to a smaller correction factor. Thus 
we need to better understand the three dimensional distribution of Abell 1689 to quantify accurately the 
underestimation of the x-ray analysis coming from a merger along the line of sight. Also it is not clear yet 
why lensing analyses need an extra mass component in the north-east (larger than the mass associated with 
the individual galaxies), whereas no x-ray emission is detected from this region. 

The lensing study of Abell 1689 is not fully satisfactory yet. On the weak lensing part, significant 
discrepancies between independent measurements of the shear profile do remain (Section 7.3) and concen- 
tration parameters larger than 20 have been claimed, which if confirmed independently will have important 
and interesting implications for cosmological models. On the strong lensing part, the very satisfactory point 
is that independent studies that used different methods and different catalogs of multiply imaged systems 
find very good agreement for the inner mass profile (Fig. [6]). However, parametric strong lensing studies 
(H06 and this work) do not give a satisfactory description of the needed second dark matter clump in t he 
north-east (whose presence has been recently confirmed bv iLeonard et al.l (|2007l ) and ISaha et al.l (J2007J)). 
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Both studies experienced some degeneracies with this second clump (see discussion in Appendix) . Therefore 
we believe that the needed mass component in that region cannot be well described with a single mass clump. 
Another evidence for this statement is that the individual RMS of the multiply imaged systems observed in 
that region (systems 6, 13, 30 and 36) are larger than the total RMS, especially for system 30. The problem 
may come from the interpretation of the strong lensing mass model. One possible explanation could be that 
this second clump is a projection effect of a large merging filament that appears in numerical simulations, 
as suggested by the optical spectroscopy of the bright galaxies in the ACS field. Clearly we need to put some 
efforts into probing the three dimensional structure of Abell 1689, in particular a lot could be learned from 
an extensive study of the line of sight velocity dispersion of the cluster members and from sz analysis. 

In this paper, we have made the well established assumption that dark matter exists. However, modified 
theories of gravity are not ruled out by this work, and we encourage people working in that area, to make 
use of the extensive data set provided in this paper, to test their theories. 

To summarize this work, we have presented an accurate mass model for Abell 1689, constructed from 
24 spectroscopically confirmed multiply imaged systems. This makes by far the most strongly constrained 
cluster to date, in terms of the number of multiply i maged systems include d with spectroscopic identification. 
Our result s are in agreement with a p revio us study by Halkola et al.1 (|2006l ) , as well as with the non parametric 
studies by iBroadhurst et al.l (|2005al ) and lDiego et al.l (|2005f ). but are much more reliable given the data we 
use. Additional spectroscopic observations are ongoing to refine the model further. We have performed a 
parallel weak lensing study, using wide field images from the CFh12k camera and devised a new method to 
select lensed background galaxies. We find very good agreement between the strong and the weak lensing 
regimes, resolving the discrepancy found in earlier work on this cluster, in particular we do not infer any high 
value for the NFW parameter . Another important result of the weak lensing analysis is that the weak lensing 
based mass estimates for galaxy clusters are reliable, even when using one band imaging with contaminated 
catalogs. With the goal of bridging the gap between the strong and weak lensing regimes, we suggest a new 
method of diagnosing the surface mass density in the intermediate region, by using the density of detected 
strong galaxy-galaxy lensing events that we report here. A detailed description of this new method will be 
described in a following paper. 



The accurate mass model presented in this work is made publicly available, thus we provide a well 
understood gravitational telescope to the community and hope it will be used to go beyond the current 
observational facilities. 
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Appendix 

A - Comparison of our Best-Fit Model with other Strong Lensing Work 

As shown in Fig. [6j we find a very good agreement with inner mass profiles derived in previous studies 
by B05 and H06. All studies also agree on the inferred concentration parameter, C200 ~ 6. 



Comparison with Non-Parametric Methods 

Non-parametric methods are feasible when the number of constraints becomes large enough, as is the 
case with the new ACS observations of Abell 1689. A clear advantage provided by non-parametric methods 
is the flexibility, i.e. the allowed range for the mass distribution is much wider than with parametric 
met hods. Abell 1689 h as been successfully modeled using ACS data and a non parametric method by B05 
and Diego et al. ( 20051 ). However, these non-parametric works were able to reproduce accurately three sets 



of unlikely multiple images (system 12, 20 and 25 of B05). Using a parametric mass reconstruction, we 
were not able to accurately reproduce these multiply imaged systems, and the parametric strong lensing 
analysis by H06 led to the same conclusion. The large freedom allowed by non-parametric methods, could 
result in the possibility of a good reproduction of partly wrong data, without a clear possibility to predict 
misidentifications of multiple images. However, the important and encouraging point is, that the overall 
mass profiles derived from parametric and non-parametric methods are in very good agreement. 



Parameterization of the Dark Matter Distribution 

On Large Scales: 

The parameters found for the two galaxy cluster scale dark matter clumps roughly agree with the ones 
found by H06. However, the description of the second north-east mass clump is not satisfactory yet. There 
is good evidence from x-ray and lensing that the dominant mass distribution in Abell 1689 is centrally 
distributed. But some extra mass is clearly required in the north-east to reproduce the observations: all 
parametric strong lensing studies need a second clump to r eproduce the multiple images configuration, even 
the very first anal ysis by Miralda-Escude &: Babul dl995 ) whic h is based on ground based imaging, and 



recent analyses by iLeonard et al.l (|2007l ) and ISaha et alJ (|2007l ) also provide evidence for substructure in 
that region. 

During the modeling of Abell 1689, we encountered some degeneracies associated with this second clump: 
if we did not restrict the value of the velocity dispersion of this second clump to be small, as suggested by 
the absence of x-ray emission from that region, we found a solution with an equally good fit describing 
this clump with a high velocity dispersion (> 1000 kms -1 ) and a large core radius of about 300 kpc. This 
solution gives roughly the same projected mass as the solution given in Table [3] since for the PIEMD profile, 
the mass scales as <Tg/r core . This large core radius is larger than the radius over which the strong lensing 
information can be found, thus a clump with such a large radius is an unlikely solution and we preferred the 
solution given in Table [3l Given the relatively large number of constraints we have, we find it puzzling to 
still find degeneracies in the description of the mass in the nort-east. 

Looking in detail at the modeling of H06 also indicates some problems in that region. H06 propose two 
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different descriptions of the second clump, one with a non singular isothermal ellipsoid, which is in agreement 
with our piemd description, and another one with an elliptical nfw profile. Depending on the model they 
use, the position changes significantly (by about 23"). Moreover, when using the elliptical NFW profile, a 
concentration parameter smaller than 1 is inferred, making this clump very flat, and it is legitimate to ask 
oneself about the reality of such a clump, since in numerical simulations, no clumps with a concentration 
parameter smaller than 3 do form. 



Recent interesting analysis by I S aha et al.l (|2007f ) gives more insight on substructures in Abell 1689: they 
use two subsets of multiply imaged systems to get two independent mass maps that agree with each other. 
Then they substracted off the best nfw fit from these mass maps, and found strong evidence for substructures 
independently from the two different data sets, in particular an extended substructure in the north-east. 
The substructures they infer are special in the sense that they are characterized by extended features much 
larger than galaxies and more massive than the stellar content, but correlated with galaxies. They call these 
extended irregular structures meso- structures that appear to be merging or otherwise dynamically evolving 
systems. 

We find the peak of the total mass dis t ributi on to coincide with the center of the x-ray emission. They 
both coincide with the BCG. IZekser et al.l (|2006! ) reported an offset between the center of the mass model 
center and the BCG of ~ 21 k pc, and H06 also r eported a similar shift of ~ 16.5 kpc. It is interesting to 
compare to the offset found in lSmith et al.l (|2005r i for a sample of x-ray luminous clusters whose properties 
are comparable to Abell 1689 (i.e. they have a high x-ray luminosity and are located at z ~ 0.2). Considering 
the irregular clusters in their sample, they find offsets between the x-ray center and the lensing center of 
mass ranging between 10 and 120 kpc. 



On Galaxy Scales: 



A similar optimization procedure has been performed by Smith et al. ( 2005h for a sample of 6 x-ray 
luminous clusters, whose properties are comparable to Abell 1689 (i.e. they have a high x-ray luminosity 
and are located at z ~ 0.2). They also included the elliptical galaxies as small scale perturbers, using similar 
scaling relations as we use in this work. They find, for a characteristic luminosity and using a piemd profile, 
(To = 180 ± 20 kms -1 and r cut = 30 kpc, in agreement with the results presented here. 



H06 and lHalkola et al.l (|2007l ) also find the elliptical galaxies in Abell 1689 to be truncated, with r cut = 
60 kpc for <7o = 220 kms" 1 (for a quite lu minous galaxy with a F775W AB magnitude equal to ~ 16.5), 
in agreement with the results of this work. Halkola et al. ( 20071 ) investigated two different possible scaling 
relations, r cut = r* ut (jj) 1 ^ 2 as assumed here, and r cut = r* ut (xr) 1 ^ 4 , corresponding to M oc L -1 / 4 . In each 
case, they found the halos to have a small spatial extent (compared to field galaxies of equivalent luminosity), 
and they were not able to discriminate between these two scaling relations, since both provide a good fit to 
the data. 

Galaxy dark matter halos in clusters have been probed in e arl ier work using galaxy -galaxy lensing 
techniques ( Geiger fc Schneider 19991 : Nataraian et al. 1998 , 2002al fbl: Limousin et al. 2007 ). These studies 
find cluster galaxies to be significantly more compact (i.e. with r cut < 50 kpc) compared to halos around field 
galaxies of equivalent luminosity (f or which no clear edge to the mass distribution has been found so far), in 
agreement with tidal stripping (see iLimousin et al.l (|2OO70 for a review of galaxy-galaxy lensing detections). 



B - New Spectroscopically Confirmed Multiply Imaged Systems 



In this section, we provide some pictures showing the new spectroscopically confirmed multiply imaged 
systems we propose in this work. Fig. [2] shows the color image of Abell 1689, with the location of the new 
multiply imaged systems. Then the following figures shows each image composing the different systems. Size 
of each panel is 27 kpc x 27 kpc. For the others multiply imaged systems, see B05 or H06. Note that the 
splitting of system 12 from B05 in two different systems has been first proposed by H06. 



Fig. 12.— Systems 12 and 31 at z = 1.83. 
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Fig. 13. — System 32 at z — 3.0. The image 32.5 in the north west is hardly detected on the color image 
but we were able to detect the two little spots on the F775W image, very close to a cluster member galaxy. 
Note that without spectroscopy, we would have been unable to reliably find this counter image. 



Fig. 14. — System 33 at z — 4.58. Image 33.2 corresponds to image 25.2 as identified by B05. 



Fig. 15.— System 35 at z = 1.52. 



Fig. 16.— System 36 at z = 3.0. 



Fig. 17.— System 40 at z = 2.52. 



